Spin tunneling and phonon-assisted relaxation in Mni 2 -acetate 
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We present a comprehensive theory of the magnetization relaxation in a Mni2-acetate crystal in 
the high-temperature regime (T > 1 K), which is based on phonon-assisted spin tunneling induced 
by quartic magnetic anisotropy and weak transverse magnetic fields. The overall relaxation rate as 
function of the longitudinal magnetic field is calculated and shown to agree well with experimental 
data including all resonance peaks measured so far. The Lorentzian shape of the resonances, which 
we obtain via a generalized master equation that includes spin tunneling, is also in good agreement 
with recent data. We derive a general formula for the tunnel splitting energy of these resonances. 
We show that fourth-order diagonal terms in the Hamiltonian lead to satellite peaks. A derivation 
of the effective linewidth of a resonance peak is given and shown to agree well with experimental 
data. In addition, previously unknown spin-phonon coupling constants are calculated explicitly. The 
values obtained for these constants and for the sound velocity are also in good agreement with recent 
data. We show that the spin relaxation in Mni2-acetate takes place via several transition paths of 
comparable weight. These transition paths are expressed in terms of intermediate relaxation times, 
which are calculated and which can be tested experimentally. 



PACS numbers: 75.45. +j, 75.30. Gw, 75.50.Tt, 75.30.Pd 
I. INTRODUCTION 



The magnetization relaxation in the molecular magnet 
Mni2-acetate with chemical formula, [Mni2(CH3COO)i6 
(H 2 0) 4 0i 2 ]- 2CH 3 COOH • 4H 2 0, (henceforth abbrevi- 
ated as Mn.12) has .ajtJbjacted much recent interest since 
several experimentsEJtmoO have indicated unusually long 
relaxation times — about two months at a temperature 
of about 2 K|-j™as well as pronounced peaks in the re- 
laxation timeEfLm in response to a varying magnetic field 
H z when applied along the easy axis of the Mni 2 crys- 
tal. These peaks correspond to an increased relaxation 
rate of the magnetization of Mni 2 and occur when H z is 
tuned to multiples of about 0.44 T. According to earlier 
suggestionsEEj this phenomenon has been interpreted as 
a manifestation of resonant tunneling of the magnetiza- 
tion, often referred to as macroscopic quantum tunneling 
(MQT). A qualitative explanation goes as follows. From 
the microscopic point of view a Mni 2 cluster acts like 
a giant spin with length s = 10 as long as the external 
magnetic field is small compared to the exchange inter- 
actions between the Mn ions, which is fulfilled in the 
experimental range considered in this paper. The relax- 
ation rate of the magnetization increases at field values 
where the spin states become pairwise degenerate. It is 
this degeneracy that determines the resonance condition. 
As the external field H z is moved away from a resonance 
the spin states are no longer perfectly degenerate, and 
therefore the tunneling probability becomes smaller and 
thus the relaxation rate. Since the spin system couples 
to the environmental phonons of the Mni 2 crystal, the 
energy levels of the spin states are smeared out. This 
leads to homogeneously broadened resonance peaks that 
are of Lorentzian shape. There are also other sources 
which lead to broadenina-af the resonances, such as hy- 
perfine and dipolar fields Hj They give rise to inhomoge- 



neous broadening with Gaussian-shaped peaks@§ How- 
ever, this stands in contrast to the measured-resonance 
peaks, which are nearly perfect Lorentzians.tj Further-, 
more, the width of the hyperfme induced GaussianaLj'E-a 
turns out to be smaller for T > 1 K than the width of the. 
Lorentzians obtained below and seen in the experiment .cB 
Similarly, dipolar interactions ttave been ruled out by ex- 
periments on diluted samples. lj Thus, for temperatures 
T > 1 K we can safely neglect hyperfme and dipolar 
fields, and the dominant source of the peak broadening 
can be explained consistently by spin-phonon effects only. 

In a || critical ..comparison between mpjdei, calcu- 
lationJ jlJ^ 'BEN and experimental dataffiM Fried- 



man et alJEi point out that a consistent explanation of 
the experimentally observed relaxation is still missing. A 
good starting point for theoretical calculations has been 
formulated by Villain et 0,1,13 where the relaxation is de- 
scribed in terms of spin-phonon interaction and a general- 
ized Orbach process. However, this approach does not in- 
clude the dependence on the external field H z . Also, one 
of the main challenges for theory is to explain the overall 
shape of the relaxation curve as well as the nearly perfect 
Lorentzian shape of the measured resonance peaksO 

In this work we perform a model calculation of the 
magnetization relaxation which is based on phonon- 
assisted tunneling. We present a self-consistent theory 
which is for the first time in reasonably good agreement 
both with the overall relaxation rate ^including all res- 
onances) measured by Thomas et alxl (see Fig. ||) and 
with the Lorentzian shape of the first resonance-peaks 
(see Figs. [?] and ^) measured by Friedman et alE3 with 
high precision for four different temperatures. 

Our model, which is introduced in Sec. ||, contains 
five independent parameters: three anisotropy constants 
A ^> B 3> B4, the misalignment angle 6 (angle between 
field direction and easy axis, the latter being taken along 
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the z axis) , and the sound velocity c. The anisotropy con- 
stant -B4 and the angle are responsible for the spin tun- 
neling. This will be explained in Sec. [II. Moreover, we 
derive the spin-phonon coupling constants in Sec. ||. It 
turns out from our calculations that these constants can 
be expressed in terms of the anisotropy. JL The constants 
A,B,B4 have already been measurecH'E2l and are known 
within some experimental uncertainty. We achieve opti- 
mal agreement between our theory and data if we pro- 
ceed as follows. In accordance with Ref. [ll] we set 9 = 1°, 
while the constants A,B,B/± are fitted to the relaxation 
data by observing, however, the constraints that A, B, B4 
are allowed to vary only within the range of their experi- 
mental uncertainties. The sound velocity c has not been 
directly measured yet (to, our knowledge). However, spe- 
cific heat measurements^ yield the Debye temperature of 
Mni2 , from which a sound velocity can be deduced that is 
in excellent agreement with our fit of the sound velocity 
c = (1.45-2.0) x l^WajWSec. |y|). Thus, in contrast 
to previous resultsliSliJO'Hlij our theory is in reasonably 
good agreement not only with the relaxation dataElO but 
also with all experimental parameter values known so far 
(see Figs. ||, [?], and^J). In addition, new predictions are 
made which can be tested experimentally: the sound ve- 
locity c and the intermediate relaxation times r n , as well 
as satellite peaks. 

extending previous work,l 12 H 14 l'E-]fc' 



In Sec. 
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U'l ti we 

make use of a generalized master equation which treats 
phonon-induced spin transitions between nearest and 
next-nearest energy levels as well as resonant tunneling 
due to quartic anisotropies and transverse fields on the 
same footing, which results in the Lorentzian shape of 
the resonances. We derive the effective linewidth of the 



Lorentzian peaks (see Sec. VI) as well as a generalized 
formula of the tunnel splitting energy (see Sec. |l|). In 
Sec. IV, we obtain the relaxation time byexactly diag- 
onalizing the master equation. In Sec. |V|, solving the 
master equation analytically, we identify the dominant 
transition paths (see Figs. |^ and [n|) and show that the 
magnetization reversal is not dominated by just one sin- 
gle path but rather by several paths which can be of com- 
parable weight. We finally note that some of the results 
of the present paper have been published in Ref. |2^ in a 
short and less general form. Here we present details of the 
derivation of these results and generalize them in various 
ways, leading to new results such as satellite peaks in the 
overall relaxation curve, relaxation time of an individual 
relaxation path, an analytical expression for the effective 
lincwidths, and a generalized tunnel splitting formula. 



II. MODEL 

In accordance with earlier work0H0lll'@0 we use 
a single-spin Hamiltonian H = 7i& + TLz + H sp + "Ht 
including spin-phonon coupling. This model turns out to 
be sufficient to describe the behavior of the Mni2-acetate 



molecule (for temperatures T > 1 K). In particular, 

H a = -ASl - BS 4 Z (1) 

represents the magnetic anisotropy where A ^> B > 0. 
The anisotropy — ASl IS depicted in Fig. |l|. We define 
the easy axis to lie along the z direction. 



A Energy 




FIG. 1. Anisotropy energy —Am 2 — Bm 4 . 

Here, S is the-,spin operator with s = 10, and A/Ub = 
0.52 -pJ0-56 fcP'B and B/k B = (1.1 - 1.3) x 10~ 3 K 
(RefsEJ andEj) are the anisotropy constants (fee is the 
Boltzmann factor). The Zeeman term 



"Hz = g^BH z S z 



(2) 



describes the coupling between the external magnetic 
field -Ha— and the spin S. The g factor is known to be 
g = 1.9B 

We denote by \m), —s < m < s, the eigenstates of 
W a + T~tz with eigenvalue 



e m = —Am — Bm + g/iBH z m. 



(3) 



If the external magnetic field H z is increased, one ob- 
tains doubly degenerate spin states whenever a level m 
coincides with a level m! on the opposite side of the well 
(separated by the barrier given by A). The resonance 
condition for double degeneracy, i.e., e m = e m <, leads to 
the resonance field 



h: 



[A + B (r 



,'2 



)] 



(4) 



As usual, we refer to n = m + m 1 — even (odd) as even 
(odd) resonances. 
The Hamiltonian 



1 



B A (Sl + St)+gfx B H x S x , 



(5) 



makes tunneling between S z states possible, where S± — 
S x ± iS y , and B^ is the fourth-order anisotropy con- 
stant. H x = IHjsin^ is the transverse field, with 8 
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being the misalignment angle. H x is assumed to be 
much smaller than H z , ri^., 9 is at most a few de- 
grees. From experiments^!! it is known that B^/hs — 
(4.3 — 14.4) x lCCi 5 K. Finally, the most general spin- 
phonon couplingc3 which is allowed by the tetragonal 
symmetry of the Mni2 crystal in leading order is given 
by 



1 



1 !)!) 



)®{Sl-Sl) + -g 2 e xy ®{S x ,S y } 



®{S X ,S Z } + 6y Z ®{Sy,S Z }) 

+ -gi{u xz ® {S x , S z } + u yz ® {S y , S z }) 



(6) 



we find (to leading order in ui a p) that the easy axis term, 
—AS 2 , is transformed into 



A(uj xz {S x , S z } + u yz {S v , S z }). 



(13) 

then yields 



Comparison with the last term in Eq. 
54 = 2A 

If the rotation matrices R a , a = x, y, z, are expanded 
up to second order, one finds terms that include symmet- 
ric elements of the strain tensor e, 



Rx 



1 

l-\8<& 







(14) 



:ffi(exa 



tyy) ® (S 2 + + S 2 _) + -g 2 e xy ® {Si - S\) 



1 



93 [(« 



{S+, S z } 



+ (e xz + iey Z ) (g) {S-,S Z }] 

+ -gi[{u xz - iuj yz ) <8> {S+, 5 Z } 

+ k z + iu) yz ) ® {5-, S z }] , 



(7) 



where i — 1,2, 3, 4, are the spin-phonon coupling con- 
stants, which we shall determine in the following. 

The linear strain tensor is defined by e — Vu, where 
u(a;, y, z) is the displacement field. Symmetrization of 
the strain tensor yields 



1 / du a dufA 



£q/3 " 2 V df3 ' da J' 

while the antisymmctrized linear strain tensor reads 

' du a 
d(3 



da 



(8) 



(9) 



with a, (3 = x, y, z. To determine gi occurring in Eq. 
we follow Dohm and Fulde.E3 The displacement 

u = 8<j) x x 



(10) 



(in leading order) is generated by rotation only. The 
infinitesimal rotation angle can be calculated by acting 
with V x (with respect to the position x) on both sides 
of Eq. (jTT 



=ivxu= 



U)y z 
U zx 
^xy 



(11) 



Applying infinitesimal rotations on the spin vector S 
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s z 



Sy U) xz U)y Z S X + Uly Z S Z 

^xzSx WyzSy S z 



(12) 



Ry = 



R, = 
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1 
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1 



(15) 



(16) 



Now we obtain from u = R z R y R x x. — x 



u = 6<p x x 



X 

Hi) y 



(17) 



By keeping derivatives of 84> a , up to second order we find 
8(b x = e xx — e yy — e zz , and cyclic permutation of (x,y,z). 
After inserting the rotated spin vector R x R y S into 
-AS 2 = -A(S 2 - Si - Si) we get for the right-hand 
side 



A(e x 



1 ii y 



)(S 2 x -S 2 )+0(e 2 ), 



(18) 



where we retain only terms that induce spin transitions. 
Comparing with the spin-phonon Hamiltonian (|TJ) one 
sees immediately that g\ — A, and thus 



gi = g 4 /2 = A. 



(19) 



Thus the coupling constants g± and g 4 are explicitly ex- 
pressed in terms of the anisotropy A. 

Finally, we note that the terms in Eq. (JgJ) that are pro- 
portional to gi^2 produce second-order transitions with 
Am = ±2, while the ones proportional to pro- 
duce first-order transitions with Am = ±1. Thus, 
Eq. ([l9]) implies that first-order and second-order tran- 
sitions are equally important fpe the relaxation. In fol- 
lowing Abragam and Bleaney,E3 it is now very plausi- 
ble to adopt the approximations \g 2 \ ~ g\ = A and 
|*73 1 Ri g 4 — 2 A (the sign is irrelevant for the transition 
rates calculated below). 
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III. MASTER EQUATION INCLUDING SPIN 
TUNNELING 

A. Generalized master equation 

In this section we derive a master equation that de- 
scribes the relaxation of the spin due to phonon-assisted 
transitions including resonances due to rfcupneling. For 
this we make use of a standard formalismE3l23 suitable to 
describe a system (spin) coupled to a heat bath reservoir 
(phonons), the latter of which is in thermodynamic equi- 
librium described by the canonical density matrix p p h for 
free phonons. That means we start from the full Hamil- 
tonian H = Ho + H p h + H sp , where Ho = H a + Hi + Ht 
represents the system, H p h the phonon heat bath, and 
H sp given in Eq. (||) is of the form H sp = J2 t Qi ® F, 
where Qi is a spin operator and is a phonon operator. 
The generalized master equation in the interaction pic- 
ture (J) in Born and Markoff approximation reads [see 
Eq. (8.1.22) in Ref. M 



dt" 



x{[Qi(i), <&(*-*")/(*)] W)Fi) 

- [Q<(*),/(t)Qi(*-*")](W(*"))}i (20) 

where Q l {t) = e mot Qie- iHot . Equation @ is valid 
for the situations where the correlation time r c in the 
heat bath is much smaller than the relaxation time r 
of the spin system. Indeed, the assumption is satisfied 
here since a rough estimate for thermal phonons yields 
t c ~ h/k B T ~ 1CT 11 s, at T = 1 K, whereas it will turn 
out that t > 1 s (see below). In this case, the integral 
kernel gives a vanishing contribution for times t" larger 
than t c , and thus one can extend the upper limit of the 
time integral to infinity and replace p(t — t") by p(t) (see 
also Ref. g Chap. 13). 

As our undamped Hamiltonian Ho has also non- 
diagonal elements in the \m) -basis it proves con- 
venient to formulate the generalized master equa- 
tion in the Schrodinger picture, i.e., with p (t) — 



Pit) = rW).«ol 



E 



dt" 



x{[Qi,Qj(-t") P (t)} mOFj) 

-iQi,p(t)Qj(-t")] {F]Fi{t"))} . 



(21) 



As the tunnel splitting generated by Ht is smaller than 
the level spacing of Ho, i.e., E mm > < |e m — e m >\ (see be- 
low), we can approximate the free propagator e ~( l / h ) n ot 
within the integral kernel by e~( l /' i ^ Wa+7 ^ z ' )t in the 
rest of our calculations!^ Next, we take the matrix ele- 
ments of Eq. (^) using p mm > = (m \p\ to'), p v 



bin 



— ~-"0 f Ifhlfh \**"J/*| v /5 r lib 1—1 I f'l l 

, ' t , and with the definitionsc3'E3 



1 f°° 

Knkln = ft2 £ <™ \Qi\ fc ) ( l \Qi I ») / dt " U ^ > 

i,3 

1 f°° 



— £ C^m'fcfcm' + ^ 



mn ^ nmmn 1 nramn 



\ - r+ - r~ 

nikkm } mmm'm' mmm'm' "> 

(22) 



it follows from Eq. (Ell) that 



Pmm' — ^_ ^o] mm ' $mm' ^ ^ Pn^mn OVnm' Pmm' ? 

(23) 

where we have considered only the secular terms [i.e., the 
"coarse-seamed" derivative was taken with respect to t in 



Eq. (EM and set [,] . = <m |[,]| m' 



The relaxation 
of the magnetization is entirely based on Eq. (^3|) . 

The difference to the usual master equation is that 
Eq. (^3|) takes also off-diagonal elements of the density 
matrix p(t) into account. This is essential to describe 
tunneling of the magnetization, which is caused by the 
overlap of the S z states. 

The diagonal elements (to = to') of Eq. ( p3f ) yield the 
master equation 



Pv 



n 



[p,T<o 



£ 



PnW mn - p m ^2 W nm- (24) 



The equation for the off-diagonal elements (to ^ to'), 



Pn 



= -^ [p,no]r 



1 'mm' Pmm' j 



(25) 



can be simplified in the following way. Accord- 
ing to Eq. (0), Qi is an element of the set 
{S% , S^L,S+S Z , S Z S+, S^S Z ,S Z S^}. Hence, we see that 



mmm m 



mmm m 



0, and we get 



2 ^ 



{Wnm- + W nm ) 



-{W m +W m >): 



(26) 



where we use the abbreviation W m = ^2 n W nm . 

Evaluation of Eq. ( |22] ) leads immediately to Fermi's 
golden rule for transition rates in first quantization [see 
Eq. (8.2.3) in Ref. ||: 

Wmn = t E I ( miV i^pi nN '^ I 2 ( N ' ^) 

JVJV' 

x 5 (E N , -E N - e m + e n ) . (27) 
Explicit evaluation yields (see Appendix A), 
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A 2 s ±1 



m ± 1 , m 



) 3 

1/ 



m ± 2 , m 



17A 2 s ±2 (e 



m±2 



-1 ' 



1927rpc 5 ft 4 e^ e " 



1 



(28) 



where s±i = (s =p m)(s ± m + l)(2m ± l) 2 , and s±2 = 
(s = Fm)(s±m + l)(s = Fm— l)(s±m + 2). Tip mass density 
p for Mni2 is given by 1.83 x 10 3 kg/m 3 !! 3 ] Here, c is the 
sound velocity of the Mn crystal, which is the only free 
parameter in our theory. As already mentioned, we are 
not aware of direct measurements of c (but see below). 
Note that the transition rates W m ±i,m, W m ±2, m are very 
sensitive to variations of the sound velocity c, as the latter 
enters with the fifth power. 



B. Spin tunneling 

We include now the spin tunneling in the generalized 
master equation fl23|). Let \m) and \m!) be two eigen- 
states of H a + Hz on the left and right sides of the 
barrier, respectively. \m) and \m') are degenerate when 
SH Z = H z LTn — H z vanishes. In the presence of tunnel- 
ing, induced by Ht, the two states form (anti) symmetric 
levels split by E mm > (for 5H Z =0). By using time- 
independent perturbation theory in higher order the tun- 
nel splitting energy E mm i can be evaluated,! 3 -! 



E,, 



m,m— 1 



m— l.m — 2 



~-m E-m — 2 



V- 



m' + l, — m' 



(29) 



Note that in this expression only steps with Am = ±1 
are allowed. However, for the present purpose we need 
to generalize Eq. ( j29| ) to situations where potentials 
^mi,mi+i S ® with arbitrary steps Am = m, — m i+ i 
(to > rrii > TOj_|_i > to', i = 1, . . . , N — 1) can occur. As 
we will show in Appendix B this is indeed possible by 
using resolvent techniques, and we find 



E 



N-l 

n 



rrii ,mi4_i 



■V, 



-m i+ i 



(30) 



where N is the lowest order of the degenerate perturba- 
tion theory, by means of which-Eq. (30) has been derived 
(TVth-order secular equatioro), giving a nonvanishing 
contribution to E mm > . For the potentials V mi>mi+1 we in- 
sert combinations of terms occurring in Ht- For example, 
the anisotropy B4 leads to transitions Am = ±4, while 
the misalignment H x leads to transitions Am = ±1. The 
summations in Eq. (|3^) can be thought of as summation 
over different paths in the Hilbert space connecting |m) 
with I to'). 

Continuing the evaluation of the first part of Eq. (|2J 
we project the undamped Hamiltonian Hq by P 



n—m,m 



on the two-state system {|? 



which yields the two-state Hamiltonian in the presence 
of a bias field (see Fig. ||) 



Ht = £m I to) (to| + 



|m) (to I + (m «-> to ) (31) 



= PH a P, 



'2™ £m 

with £ m = £ m + g/is<5iJ z m and the energy eigenvalues 

E T = \ U + W ± yJ(U-U') 2 +En m > ■ pro- 

vides a valid description as long as the level splitting 
remains smaller than the level spacing, i.e., 



A = 



\J \U - U') 2 + E 2 mm , « \e m(t> - e m0)±1 \. (32) 



We have checked that between two main resonances this 
condition is satisfied for the states |tot) and \m' T ) of the 
dominant paths and for each degenerate pair of states 
\m), \m') with e m ,e m > < e mT ,e m ^ (see Sec. 0). 

Next we insert the two-state Hamiltonian Ht into the 
generalized master equation (E3|), which yields 



iE„ 



(Pr, 



) - W m pm + W mn Pn 



(33) 



E 



and 



Pn 



1 £ \ t-^mm' 



{Pm-Pm'), (34) 



where £ mm ' = £ m — £ m ', and likewise for to <-> to'. Ul- 
timately, we are interested in the overall relaxation time 
r of the quantity p s — p_ s (see Sec. |v|) due to phonon- 
induced transitions. This t turns out to be much longer 
than Td = l/'Ymm'i which is the decoherence time of 
the decay of the off-diagonal elements p mm > oc e~*/ Td 
of the density matrix p. Thus, we can neglect the time- 
dependence of the off-diagonal elements, i.e., p mm > ~ 0. 
Physically this means that we deal with incoherent tun- 
neling for times t > Td.E£l Inserting then the stationary 
solution of Eq. ( p4| ) into Eq. (|3^), which leads to the 
complete master equation including resonant as well as 
nonresonant levels, 



Pm = 

where 



-W mPm + W mnPn + T™' (p m , - p m ) , (35) 



W m + W m ' 



^l m/ + h 2 (w m + w ml y 



(36) 



is the transition rate from m to to' (indweed by tunnel- 
ing) in the presence of phonon damping.^ The relaxation 
dynamics of the resonances described by 1/r™ ~ 10 -7 s 
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(see Fig. |T^) turns out to be much faster than the phonon- 
induced overall relaxation, i.e., 1/r™ <C r > Is (see 
Fig. H). Thus, our derivation based on the assumption 
r 3> Td is self-consistent since 1/r™ ~ rj. Note that 
Eq. (|3^) is now of the usual form of a master equation, 
i.e., only diagonal elements of the density matrix p{t) 
occur. For levels k ^ m,m' , Eq. (p5|) reduces to 



Pk 



-W k p k + ^2 W knPr, 



(37) 



We note that r™ has a Lorentzian shape with respect 
to the external magnetic field SH Z occurring in £ mm '. 
The H z dependence of W m around the resonances turns 
out to be much weaker (see below) and can be safely 
ignored. It is thus this that will determine the peak 
shape of the magnetization resonances (see below and 
Figs. |3|-||). Note that in Figs. ||-|| these Lorentzians are 
truncated at the center of the peak by the spin-phonon 
transition rates W m and W m ' in such a way that the 
effective linewidth (defined as the width at half of the 
height of the truncated peak) is much larger than (W m + 
W rn i )/2. This needs some further explanations, which are 
given in Sec. VI, after we have discussed the relaxation 
times. 
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\A) = Jj(|m>-|m')) 
s ) = 75 (M + !'"'» 



FIG. 2. Tunneling configuration. 



where the elements of the vector p(t) are the diagonal ele- 
ments of the density matrix p. Within the interval I ni .n 2 
[see Eq. (|65|)] delimited by the two main resonances n\ 

and ri2 only the tunneling rates Tml and Tml [Eq. (|3^) 1 
for which mi+mj = n\ and m2+m' 2 = n 2 (see Sec. |V A| ) 
are allowed to be included for self-consistency reasons: 
the tunnel splitting ( |32| ) entering T is only valid within 
this interval I nii n 2 - If w ii i = 1,2,..., 21, are the eigen- 
values of the rate matrix W, the dominant relaxation 
time of the spin system is given by 

1 



r = max < — 



Re Wi 



(39) 



The eigenvalues uii turn out to be nondegenerate with 
the smallest one being far separated by a factor of at 
least 10 4 from the remaining ones. The result is plotted 
in Fig. [I, where the overall relaxation rate r is shown as 
a function of H z at T = 1.9 K. It is important to note 
that all the resonance peaks are of Lorentzian shape. We 
note that in our model the even resonances are induced 
by the quartic anisotropy, whereas the odd resonances 
are induced by product combinations of B^S^. and H X S X 
terms [see Eq. (|o|)] . For the plot in Fig. || we set£ = 1° in 
accordance with the experimental uncertaintyjij leading 
to a maximal transversal field H x of about 350 G. 
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FIG. 3. Full line: semilogarithmic plot of calculated relax- 
ation time r as function of magnetic field H z at T = 1.9 
K. The optimal fit values (see text) are A/ks = 0.54 K, 
B/k B = 1.1 x lCr 
andc= 1.45xl0 3 
Ref. @. 



3 K, and B A /k B = 8.5 x HT 5 K, 8 — 1°, 
m/s. Dots and error bars: data taken from 



IV. RELAXATION TIME 

A. Numerical diagonalization of the master equation 

In this section we give the results of our exact eval- 
uation obtained by a numerical diagonalization of the 
master equation. For convenience we now write down 
the master equation (|35|) as a vector equation, 



(38) 



p(t) = Wp(t), 



B. Comparison with experimental data 



For comparison wc alsip^kicludc in Fig. || the data re- 



ported by Thomas et a/.L}t_j We have optimized the fit 
(as explained in the Introduction) in such a way that the 
fits of the model parameters, given by 



A/k B = 0.54 K, 
B/k B = 1.1 x 10~ 3 K, 
B 4 /k B = 8.5 x 10~ 5 K, 



(40) 

(41) 
(42) 



G 



are roughly within the reported experimental uncertain- 
ties of Refs. ^l| and |2^ (see above). The value of B4 is in 
excellent agreement with recent measurements performed 
in Ref. [hi Our fit of the sound velocity yields 



= 1.45 x 10 3 m/s. 



(43) 



There is a difference between odd and even resonances, 
i.e., the relaxation time r at an even resonance peak is 
about three times smaller than the one at a subsequent 
odd resonance peak. It should be mentioned that almost 
identical plots are obtained for 0.5° < 9 < 3°, as can be 
seen in Figs. ||-[| The present theory holds for \H X \ < 
1000 G (which is well satisfied here), otherwise the shift 
of the levels \m) due to the perturbation H X S X must be 
taken into account. For example, for the resonance n=3 
the relevant tunneling takes place between |4) and | — 1). 
The dominant second-order shift 

(1 \gfx B H x S x \ 0) (0\g» B H x S x \l) 
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FIG. 4. Full line: semilogarithmic plot of calculated relax- 
ation time r as function of magnetic field H z at T = 1.9 K. 
Here 6 — 0.5° has been chosen. Dots and error bars: data 
taken from Ref. 



(1 \gfJ, B H x S a 



T -l — £0 

I 2) (2 \g^i B H x S x 



£1 - £2 



= k B x 40 mK 

< \sx\ =k B x 2.3 K 



clearly shows that the unperturbed states {|ro)} are a 
good zeroth-order approximation. It is also important 
to know whether the second-order shifts caused by the 
perturbation 7ix are negligible compared to the tunnel 
splitting E mTJn > . Explicitly, we find 





_ A (2) 


/k B 


= 8.5 mK 


(n = 


4), 


Af) 


-aL 2 ] 


/k B 


= 13.2 mK 


(n-- 


= 3), 






/k B 


= 5.0 mK 


(n = 


2), 


A(2) 
^2 




fk B 


= 0.5 mK 


(n = 


1)> 


A (2) 


-aL 2 1 


fk B 


= (n = 


0). 
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FIG. 5. Full line: semilogarithmic plot of calculated relax- 
ation time r as function of magnetic field H z at T = 1.9 K. 
Here 6 = 2° has been chosen. Dots and error bars: data taken 
from Ref. ^ 



X[S] 



(2) 

where Am T is the second-order shift of the unperturbed 
states |mx) of the dominant paths. These renormaliza- 
tions cause a very small shift of the resonance peaks, e.g., 
the shift for n = 3 is 0.6 mT. The relevant tunnel split- 
ting energies of the odd and even resonances are about 
the same (except E^—i): 

£4.-4 ~ #3,-2 ~ #5,-3 ~ k B x 45 mK, 
E±,-i/k B w 130 mK, E 6i _ 2 /k B w 40 mK. (46) 

For comparison, Ei,-ilk B w 1 K [see also Eq. (|66|)1. In 
conclusion, the diagonal elements of the shifts of the non- 
degenerate perturbation theory are much smaller than 
the off-diagonal elements of the shifts of the degenerate 
perturbation theory (see Appendix B). Thus our assump- 
tion of quasidegeneracy is very well satisfied. 
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FIG. 6. Full line: semilogarithmic plot of calculated relax- 
ation time r as function of magnetic field H z at T = 1.9 K. 
Here 8 = 3° has been chosen. Dots and error bars: data taken 
from Ref. Q 
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We note that there are satellite peaks in Figsj3| |[ the 
origin of which will be explained below in Sec. [V| 

In Figs. and H we plot the peaks of the first resonance 
at H z = 0, which is induced only by the B4 anisotropy, 
for four different temperatures, namely T = 2.5,2.6,2.7, 
and 2.8 K. The four peaks (like all others) are of single 
Lorentzian shape as a result of the two-state transition 
rate F"^ given in Eq. (|36|). For comparison we plat in 
Figs. [7] and || the data reported by Friedman et alh3 for 
the same temperatures (no error bars, however, are given 
in Rcf. The optimal fit values are 



0.56 K, 

1.3 x 1CT 3 K, 



A/k B 
B/k B 

B 4 /k B = 14.4 x 1CT 5 K, 
c = 2.0 x 10 3 m/s. 



(47) 
(48) 
(49) 
(50) 
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FIG. 7. Full line: plot of calculated relaxation rate V — 1/r 
as function of H z for the first resonance peak at T = 2.6 K. 
The Lorentzian shape is due to T 7 ^ in Eq. (ji^) . The optimal 
fit values (see text) are A/k B = 0.56 K, B/k B = 1.3 x 1CT 3 
K, and B A /k B = 14.4 x 1(T 5 K, 9 = 1°, and c = 2.0 x 10 3 
m/s. Dots: data taken from Ref. [ll]. 
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FIG. 8. Full lines: semilogarithmic plots of calculated re- 
laxation rate F = 1/r as function of H z for the first resonance 
peak at (a) T = 2.5 K, (b) T = 2.6 K, (c) T = 2.7 K, and 
(d) T — 2.8 K. All peaks are of single Lorentzian shape. The 
optimal fit values are the same as in Fig. (jj. Dots: data taken 
from Ref. [n]. 



Note that these values are the same for all four tem- 
peratures, which means that our peaks fit also the tem- 
perature dependence of the relaxation time. The fitting 
parameters turn out to be somewhat larger than the ones 
used in Fig. || [see Eqs. (|40|)-(^2])|, which could be caused 
by sample differences, e.g., in volume-to-surface ratio 
and/or in shape anisotropy of the samples, etc. Indeed, 
the sample of Ref. [ll] consists of many small crystallites 
in contrast to the single crystal used by Thomas et alu In 
any case, the differences are small, and the sound veloc- 
ity c seems to be within the expected order of magnitude. 
Clearly, it would be highly desirable to check this predic- 
tion by an independent and direct measurement of c. On 
the other handjiil we can get an independent estimate for 
c from the specific heat and the Debye temperature O u 
which was recently measured in Mni2.c3 The reported 
value iSr-©_D = (38 ± 4) K, and making use of the Debye 
relatione 3 ! 



k B Q B = Tujj b — hcku, with n 



I- 3 
6tt 2 



1 



we find 



c = (1.77 - 2.18) x 10 3 m/s, 



(51) 



(52) 



where u>d is the Debye frequency, ko the Debye wave 
vector, and Vq — 3716 A 3 the unit-cell volume. Com- 
parin g th is value for c with the one obtained before, see 
Eqs. ( J43] ) and (|50|), we see that the agreement is very 
good. This result corroborates not only our prediction of 
c but also our values obtained for the spin-phonon cou- 
pling constants g,. 

Finally we also mention that the prefactor 
A 2 s±i/12Trpc 5 h 4 of our spin-phonon rates [Eq. pSj)] is 
in excellent agreement with thepVjalue of the parameter 
denoted by Cria a recent paper M Note that their fit of 
the parameters C is not as precise as ours, because C is 
assumed to be independent of the spin states 

To summarize our results obtained so far, we see that 
the agreement between theory and experiment is satis- 
factory; in particular we emphasize that there is no free 
fit parameter. Thus, our model and its evaluation seems 
to contain the essential physics responsible for the mag- 
netization relaxation in Mni2. 



C. Comparison with previous results 

In comparison to previous results we obtain much bet- 
ter agreement between theory and experiment for the fol- 
lowing reasons. For this comparison we can restrict our- 
selves to the work of Fort et al.p3 since — as far as we 
are aware of — it has produced the best agreement with 
the relaxation dataEI thus far. First, the spin-phonon cou- 
pling constants g±, 52, and (73 are explicitly given-in our 
work for the first time (34 has been found beforetS). As 
shown in Sec. |l| we find them to be of order A = 0.56 K, 
and it is this value which leads to good agreement with 
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all known experimental data which invohie these coupling 
constants.ErQa In contrast, Fort et alS3 set arbitrarily 
the f/j's to values of 15 K and 30 K, which is clearly in 
contradiction to our microscopic values. Moreover, our 
value for the fit parameter Jl^ulfilk the constraints pf 
independent measurements J13't3'Ej£3 while Fort et alxB 
obtain a B4 value which is about 30 times smaller than 
the measured value in Ref. ||l]. From Figs, fj] and || we see 
that the temperature dependence of the relaxation time 
agpees quite well with the measurements of Friedman et 
al.lA. Such a fit has not been attempted before. 



V. RELAXATION PATHS 

A. Analytical result 

In order to get a better physical understanding for the 
relaxation process of the spin system it is instructive to 
determine the dominant transition paths via which the 
spin can relax into its ground state. For this we derive an 
approximate analytic expression for the relaxation time 
denoted by r* (to distinguish it from the exact r obtained 
in the previous section). First, we solve the master equa- 
tion for one particular transition path n which does not 
intersect with other paths. For H z > we find (deriva- 
tion is given below) 



1 



0/ 3( £ _ 



E 

{m,}„ 



m i+1 



(53) 



where fC- +1 = £l m ,_ mi+1 = W mi+im or T%* +1 , depend- 
ing on the particular path n characterized by the sum 
over the levels m (see Figs. || and [To]). Equation ( p33| ) 
holds for arbitrary initial (e^) and final (e/) energies, and 
for arbitrary steps Am = rrii+i — m$ (see below). 

We now turn to the derivation of the relaxation time 
of a cascade including the external field H z . For this we 
need to go beyond the results obtained previouslyt3 for 
H z — 0, which requires a non-trivial extension. We start 
with rewriting the master equation (Bq) as 



Ps — ^to-lAth 

,, j+1 Pm j ^ 



iL mj Pmj "m, Cm,: 



P- 



(54) 



with rrij e ]— s, s[, rrij > rrij+i, j = 1, . . . ,p < 2s — 1, 
and p mj — dpjnj/dt. We consider now the stationary 
limit of Eq. (p4) which we define by 



Ps = - J, p mj = 0, p-s = J, 



(55) 



where the first and last equation express conservation of 
the probability current J, which we assume positive for 
H z > and independent of to. Equation (^5|) leads to 
p + 1 equations, 



T(t) — fl" H+1 n - fl mi n 
uyu)—!,L mi p mi " mi+ irm i+ u 

and by solving for p mi+1 we get 



Q mi J 

_ "m i+ i J 
Prrii — --.mj+i Pm i + 1 + ^.m i+1 



(56) 



(57) 



where we have introduced i = 0, . . . ,p < 2s— 1, and too = 
s. To simplify the following treatment we assume detailed 
balance also for the tunneling processes. This approxi- 
mation has little effect on the final result which turns out 
to agree very well with the exact relaxation time r where 
no such approximations are invoked. Inserting then the 
detailed balance relation = e /3Cem *+i _Cm « ) 

one obtains 



Prrii — c Prn i+1 i „m, + i 



(58) 



In order to get an equation that depends on p_ s and p s 
only one has to sum over the following p + 1 equations: 

Ps — e p mi -r , 

e 0(e mi -e s ) _ e f3(e m2 -e s ) , J J_ \ 

e Pm-t — c pmi \ ' 

JHsmt-e.) - e /3(e m3 - £s ) Je 0(e m2 -e 3 ) 

& fJm 2 — e f J rn 3 \ O m ' 3 ' 



p P(e mp -e- s ) _ D{e. 



a 'mi 



(59) 



In the special case of H z = QtJ.c., e s — e_ s = 0, Eq. (|5£ 
agrees with previous results. ta 

Taking the time derivative of Eq. (|59|) and using 
P-s/Ps = — 1 we find 

j e /3(e m ,--£a) 



and thus 



mi 1 

= -2J. (61) 
The solution of the last differential equation is 

J(t) = J e- t/T \ (62) 
with the relaxation time 

T* = a, m-^ , H z > 0. (63) 
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Finally, the summation ^{mjk m Eq. ( |53| ) is denned as 
the summation $^ m . in Eq. (|63| ) taken only from m; n itiai 
to m final + 1, where |m ini tiai) , |mfi na i) denote any two 
neighboring vertices (where paths intersect) in Figs, [j], 
|§ [[§, and |o] (see below). 




|10) 



|-io) 



FIG. 9. Spin relaxation paths (from m = 10 to m = —10) 
for < g[iBH z < A + 13B. Full lines: thermal transitions 
due to phonons. Dashed lines: dominant tunneling transitions 
due to B4 and H x terms. Dotted lines: tunneling transitions 
that lead to satellite peaks [included in the numerical diag- 
onalization of the master equation (jig)]. The states where 
paths intersect are denoted as vertices. 




FIG. 10. Serially reduced diagram associated with Fig. |S)| 
In order to understand the analytical evaluation of the re- 
laxation diagram in Fig. |^ better, tunneling transitions that 
lead to satellite peaks are excluded. The relaxation times 
t„ are given in Eq. (||). For \H X \ < 0.05 T only the path 
n — > T3 — *■ T7 is dominant (see Figs. ^ and |^). 



Similarly one can solve the rate equations (55) for J < 
H z < 0. Then we obtain 



1 



1 -|- e P(e s -c- s ) 



s /3(e mi -e_ a ) 



H z < 0, (64) 



which for H z — (i.e., e_ s = e s ) and steps Am = ±1 
reduces to the result found in Ref. |lj. 

If there is more than one path contributing to the re- 
laxation (which is typically the case in the region between 



two resonances), we have to account for intersections at 
vertices. For this we associate with each path a proba- 
bility current J n = p n , and interpret Eq. (53) in terms 
of a serial circuit with the summands playing the role of 
"resistances." This allows us then to set up flow diagrams 
for J„ (see Figs. ||, [l(], and[l£ -21), which obey the analog 
of Kirchhoff's rules: 

(Kl) J n — 0: The sum over all incoming and out- 
going currents vanishes at a vertex (current conserva- 
tion) . 

(K2) J2 n JnTn = AN: The sum over all voltage drops 
{JnT n ) is equal to the source-drain voltage AiV = p s — p- s 
for any closed path (probability conservation). 

The total probability current is given by J = AN. For 
every interval 



■Mil ,«2 



mT,im T 1 



HT m ' 2 \ (65) 

[see Eq. (Q)], where n% = tot.i +m' T 1 , n 2 = tot. 2 +m' T 2 , 
< m — n 2 + 1 < 3, a set of equations is given by 
the rules (Kl) and (K2). For every set we derive the 
relaxation time r* „ = AN/ J. 

Figure p[ shows the complete, Fig. IICI its serially re- 
duced flow diagram for < H z < -±- (A + IW). From 
(Kl) we get 



J = 


Jl- 


h J2, 


J2- 


-h 


— Je , 


Jl = 


J3- 


h J A , 




-h 


= Jr, 


Ji = 


J5- 




J7- 


- Js 


= J, 



while from (K2) we get 

AN = Jm + J3T3 
J3T3 = J4T4 + J5T5 
J2T2 = Jm + J4T4 

JsTs = J5T5 + J§T§ 

From these equations we obtain 
7-0,1 i H z 

T4T7T5T2 " 

T S T 7 T 2 T 3 - 

T~8T 3 T 2 T§ - 

T8T3T5T6 - 

T8T 3 riT 6 - 

T7T1T5T6 - 

T8T3T1T2 - 

T 8 T 7 T 2 T e - 



J7T7 , 

J5T5 , 
J7T7 . 



= (T4 7i r 5 r 2 + r 8 t 4 n r 2 

T4T 7 T 2 T 6 + T S T 7 T 5 T 6 + 



T4T 7 T 2 T 3 + T$T4T 3 T 7 
T8T4T3T2 + T~4T 3 T 2 T 6 
T%T 7 T^T2 + T~9,TiT 2 Tq 
T~7TlT 5 T 2 + T 3 TlT 2 T 6 
T 7 TlT 2 T 4 + T S T 6 TlT 7 
T7T1T2T3 + T S T 3 TlT 7 
T 8 T4T 7 T 6 + T4T 3 T 5 T 6 
T4TlT 5 T 6 ) J (t S T 5 T 2 + T 8 T 5 T 6 + T 8 T 3 T 2 + T 8 T 4 T 3 -\ 
T8T4T e +T 8 T4T 2 +T4T 2 T 6 +T4T 7 T 6 +T S T 3 T 5 +T 7 T 2 T 6 
T 7 T 5 T 2 +T 7 T 2 T 3 +T4T 1 T 5 +T 3 T 1 T 5 +T S T 1 T 5 +T 8 T 3 T 1 



T$T 7 T 2 T4 -+ 

T 8 r 3 r 5 T 7 4 
T 8 r 4 r 3 T 6 4 

T S TlT 5 T 6 4 
T 3 TlT 5 T 7 4 

T4T1T5T7 4 

T 3 TlT 5 T 2 4 
T4T3T5T2 4 



T S T4TlT 6 - 
T4T7T 5 T 6 - 
T4T3T5T7 - 
7"8T'3T'5T2 - 
T3T1T5T6 - 
TYT1T2T6 - 
T 8 T4riT7 - 

TST1T5T7 ~ 
T4TiT 2 Tq - 

F r 4 T 5 T 6 - 
' t 7 t 5 t 6 - 

T 3 T 2 T S - 



T 3 T 5 T 6 + T 3 T 5 T 2 + T 3 T 5 T 7 + T4 T 3 T 7 + T4 T 3 T 5 + T 8 Tq T\ + T 3 T\ T 7 + 

T7 Ti t 6 + T4 n t 6 + t 3 n t 6 + T4 n t 7 4- t 8 r 4 n + r 8 r 2 t 6 + r 4 r 5 t 2 + 
T4T 3 r 6 4- r 2 r 4 T7 4- t 7 tit 5 ). 

When Tq -4 is plotted as function of H z there is no vis- 
ible difference between the exact t obtained in Sec. IV 
and this approximate t*, which confirms that the dia- 
gram in Fig. ^| contains the physically relevant relaxation 
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paths for the interval Io,i- Similar results are obtained 
for the other intervals, whose diagrams and calculations 
are shown in Appendix C. 

Finally, near a resonance [|<5iJj < w', see Eq. (74)] 
the above expression for r* [Eq. (p3[)] strongly simplifies 
since we find that there is only one dominant relaxation 
path which involves only one tunneling channel. This fi- 
nally explains why the peak shape is given by a single 
Lorentzian. We call the five strongest broadened reso- 
nances in Figs. |^-|^ the main resonances. For every main 
resonance n we have identified [using Eq. (j53|)] its dom- 
inant path and its associated tunneling channel between 



im 2 m 2 
mimi' 



the states |tot) and 



These states are 



rriT m\ 



(66) 



Our calculation of the intermediate relaxation times r ra 
provides a further prediction which could be tested with 
NMR techniques of the type described in Ref. H 



B. Satellite peaks 

Beside the main resonances there are also other nar- 
rower resonances (see Figs. |[^) that are a direct con- 
sequence of the fourth-order anisotropy constant B [see 
Eq. (H)]. Indeed, if the plots around one peak are magni- 
fied further, satellite peaks become visible (see Figs. |lT| - 
[Ti] ) . In order to understand the occurrence of these satel- 
lite peaks it is instructive to look at Fig. [l8| below. There 
are several paths which can be used in the relaxation 
process. As we include the fourth-order anisotropy term, 
—BSz, the resonance condition is not the same for every 
level [see Eq. (|j)]. Hence, very narrow peaks show up, 
which can be seen only at high resolution. In Fig. |l8| 
several additional tunneling paths, some of which are re- 
sponsible for the satellite peaks in Figs. [□] and [l^, have 
to be drawn (represented by the dotted lines in Fig. 18). 
For example, the tunnel splitting energy of the path from 
|4) to |— 2) is proportional to H X B±H X (third-order per- 
turbation) , where the ordering of the factors corresponds 
to the chosen path. Due to the presence of H x the width 
of the satellite peak (see next section) depends on the 
misalignment angle 9. If one takes a close look at our 
high resolution plots this difference between Fig. |ll] and 
Fig. [lj is observable. It must be noted that we consider 
only tunnel splitting energies up to second order in B4 
and third order in H x (also combinations such as B\H^) 
for all the main and satellite peaks. Narrower satellite 
peaks are neglected.c3 The distance between a 

satellite peak and its associated main peak caused by a 
main resonance is given by Eq. (kf), 



rrmimi 

nB 



g^B 



HI 



2 , it 2 /2 

1 mi — m 2 — m2 



(67) 



where mi, mi' {m 2 , m 2 ) are responsible for the satellite 
(main) peak, and n — mi + mi' — m 2 + m 2 . It would 
be interesting to search experimentally for these satellite 
peaks, which requires a higher experimental resolution of 
the peaks than achieved so far. 

X [S] 




H Z [T] 



FIG. 11. Full line: semilogarithmic plot of calculated re- 
laxation time r as function of magnetic field H z at T = 1.9 
K in the interval 3A/2g/j,B < H z < 5A/2g/j,B with a higher 
resolution. The tunneling transition from |5) and |— 3) is re- 
sponsible for the main peak. Two satellite peaks are visible. 
The left (right) one is due to the tunneling channel between 
1 4) and |-2) (|6) and |-4)). Here 6 = 2° has been chosen. 
Dots and error bars: data taken from Ref. 



x [s] 
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H Z [T] 



FIG. 12. Same plot as in Fig. [d], but with a misalignment 
angle of 9 = 3°. Dots and error bars: data taken from Ref. [?|. 
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FIG. 13. Full line: semilogarithmic plot of calculated re- 
laxation time t as function of magnetic field H z at T = 1.9 
K in the interval 7A/2g/j,B < H z < 9A/2g/j,B with a higher 
resolution. The tunneling transition from |6) and | — 2) is re- 
sponsible for the main peak. Two satellite peaks are visible. 
The left (right) one is due to the tunneling channel between 
|5) and | — 1} (|7) and |— 3)). Here 6 = 2° has been chosen. 
Dots and error bars: data taken from Ref. |^. 



x [s] 



10 4 
10 3 
10 2 
10 




H Z [T] 



1.5 1.6 1.7 1.8 1.9 

FIG. 14. Same plot as in Fig. |l3| but with a misalignment 
angle of 8 = 3°. Dots and error bars: data taken from Ref. [t| 



VI. WIDTH OF THE LORENTZIANS 

In this section we give a physical interpretation of 
the effective half-width of the Lorentzian peaks in our 
plots. In order to get an expression for the width of 
our main and satellite peaks consider a Lorentzian T{H Z ) 
with linewidth w (see Fig. |l5| ). If the upper part of this 
Lorentzian is cut off (where the curve is already very nar- 
row) and both ends are connected by a horizontal line one 
obtains a curve that still has the same single Lorentzian 
shape for all practical purposes but now with an effective 
linewidth w' > w. Changing the tunnel matrix element 
E m m' results in a different truncation of the Lorentzian, 
thus changing the effective linewidth w' . We shall now 
estimate the effective linewidth w' and compare it with 
the one obtained from the exact 1/r. Taking only the 
largest terms of Eq. (m) gives a rough approximation of 



the relaxation time near a resonance where the states |m) 
and |m') are degenerate, 



1 



1 + e0(*-.-*.) 

e (3(s m -e„ s ) 



3 (3{e m + 2-e a ) 



/9(e ro /-e.) 



m,m+2 

-) 

pm' / ' 

m / 

Using the detailed balance relation 

Wm,m+2 _ e /3(e m+2 -e m ) 
Wm+2 .in 

we obtain the following approximation 

e /3(e m+ 2-e s ) / 2 



w, 



m' — 2,m' 



v m,m+2 



1 



(68) 



(69) 



(70) 



where we assumed that W m ,m+2 ~ W m >-2,m'^ In the 
limit £ mm ' — * the phonon-damped tunneling rate F™ 
is much larger than W m ,m+2, so 



lim t 



2 e /3(e m+ 2-£s) 



(1 



w, 



(71) 



m,m-t-2 



The half- width of t'(H z ), denoted by w' , is then deter- 
mined by the condition t'(w' /2) = t'(0)/2. This condi- 
tion is fulfilled when 



pm 



w, 



m , m + 2 



(72) 



Thus we obtain the expression for the effective linewidth 



\m - m'\gn B 



El 



h 2 (Wm+W m >) 



-,1/2 



ir, 



m,m+2 



(73) 



Since the height E 2 nm ,/h 2 (W m + W m >) of the Lorentzian 
L™ is very large compared to its linewidth |m — 
m'\gnBw/h = (W m + W m >)/2 and W m + W m - « 
2W m ,m+2 for the dominant paths (see Sec. 0) we get 
the following reasonably accurate approximation for the 
effective lincwidths in our plots: 



2 3/2p; 



\m — m f \gfiB 



(74) 



Comparison with our exact calculations of the relaxation 
time shows that w' of Eq. ( f74|) gives a very good estimate 
for the effective linewidth of the peaks in our plots (see 
Figs. §-| and HHH). 
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FIG. 15. Truncated Lorentzian TJ 4 with A/k B = 0.56 
K, B/k B = 1.3 x 10" 3 K, and B 4 /k B = 14.4 x 10" 5 K, 
6 = 1°, and c = 2.0 x 10 3 m/s, w' = 37.4 mT (w' agrees 
very well with the width of the Lorentzian in Fig. Q), and 
Fj 4 (-u//2) = 2.4 x 10 6 i 
the dashed line. 



The truncation is indicated by 



VII. CONCLUSION 

We have presented a comprehensive theoretical de- 
scription of spin relaxation due to phonon-induced tran- 
sitions and tunnel resonances. Deriving a generalized 
master equation (in Born and Markoff approximation) 
we obtain an exact numerical evaluation of the overall re- 
laxation time r as function of the longitudinal magnetic 
field H z comprising Lorentzian-shaped peaks. In order to 
perform this evaluation we calculate the phonon-assisted 
transition rates of the spins, the spin-phonon coupling 
constants, and the tunnel splitting energy, for which a 
generalized formula is derived. The fourth-order diago- 
nal terms in the Hamiltonian give rise to satellite peaks, 
the experimental observation of which requires a higher 
resolution of t(H z ) than achieved so far. Our approx- 
imate analytical solution of the master equation yields 
a clear physical understanding of the relaxation process 
by revealing the relaxation paths that are followed by 
the spin. This solution provides the prediction of all in- 
volved intermediate relaxation times r„, which can be 
tested experimentally. The results of our model calcula- 
tion agree well with all known data. For the first time 
we have been able to get agreement between theory and 
the entire relaxation curve. In addition, we have obtained 
reasonable agreement between theory and four single res- 
onance peaks recently measured to high accuracy at four 
different temperatures. The formalism presented in this 
work has been applied to the specific parameter values of 
Mni2, but many results derived here are generally valid 
and can be used for similar spin systems as well. 
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APPENDIX A: SPIN-PHONON RATES 

In order to evaluate the spin-phonon rates W mn of 
Eq. ( p^ ) we first change to the Fourier representation. 
If q is a phonon wave vector, we can write u(x) as fol- 
lows: 



u(x) = ^=£u(q)e^ x , 



(Al) 



with N being the number of unit cells. Hence 



s(x) 



E 



q x u x (q) q x u y {q) q x u z {q) 
q y u x {q) q y u y (q) q y u z {q) 
q z u x (q) q z u y {q) q z u z (q) 



(A2) 



After (anti)symmetrization, these matrix elements can 
be inserted into the expression (0), 



= -jL X/ l^9i[qxU x {q) - q y u y {q) 

3,1 ^ 



®{Sl + S 2 _) 



+ -g2[q x u y {q) + q y u x {q)\ ® (Si ~ S 2 + ) 

+ ^93[qxU z (q) +q z u x (q) 

~ i(q y u z (q) + q z u y (q))} ® {S+, S z } 

+ ^93[qxU z {q) +q z u x (q) 

+ i(q y u z {q) + q z Uy(q))} ® S z } 

+ ^94[qxU z (q) +q z u x (q) 

- i(q y u z (q) + q z u y (q))} ® {S+, S z } 

+ ^9i[qxU z {q) - q z u x (q) 



+ i(%uM - Wq))] ® {S-,S z }}e 1 ^. (A3) 

Rj are the positions of the Mni2 molecules. 

We proceed with the canonical transformation 
(u,p) — * (ct,c). c^) = eqCq" 1 annihilates (creates) a 
phonon with wave vector q and polarization e q , and 



where M is the mass per unit cell. Inserting Eq. (A4) 
into Eq. (A3) and considering only the spin of the Mni2 
molecule at = yields 
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n 



sp 



E 



2iVAfw 



x {~<?i[fe(4 + cx) - %(4 + c y )] ® [S% + S 2 _) 

+ £<?2[<Z*(4 + c„) + q y (ci + c x )] ® {S 2 _ - S 2 + ) 

+ ^9s[(qx~ iq v )(cl + c 
+ <7z(4 + c x ~ *4 - i^) 

+ g#3 + (4 + c 
+ ? 2 (4 + + ^4 + 

+ ^.94 [(q x - igj,)(4 + c 

- q z (cl + c x - ic\ - ic y ) 

+ ^94[(<lx +iq y ){cl +c z 

- q z {c x + c x + id. + icy) 



{S+,S Z } 
{S-,S Z } 
{S + , S z } 

{S^,S Z }}. (A5) 

This expression can be used to evaluate the transition 
probability. We employ the following standard relations: 

c |n) = \fn \n — 1) , 
c) \n) = Vn + 1 \n + 1) , 



S- |s, m) = y/ (s + m)(s — m + 1) |s, m — 1) , 
S+Isjto) = \/(s — m)(,s + m + 1) |s, m + 1) . (A6) 

The transition rate W_2 = W m -i,m [see Eq. ( p7|) 1 for 
m — > m—2 (e m -2 ^ £m) can now be calculated in second 
quantization (n a = n^ a , a — x,y,z denotes the num- 
ber of phonons with wave vector q, polarization mode A, 
and oscillation direction a, and the thermal average over 
phonons is left implicit), 



2tt 



W_ 2 = y^| (n q , =F l,m — 2 |W sp | n qV . 
q' 



|2r/ 



.9? 



4 {q x (n x T l^K/ - q y (n y T 1|4 +) K 
x \(m- 2\S 2 _\m)\ 2 



x |(m - 2\S 2 _\m)\ 

1 7TS_ 2 



NMu} 



gliix - qy) 1 



S± 

l 
o 
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(fe + 9y) 2 



(A7) 



where s_2 = (s + m)(s — m + l)(s + to — l)(s — m + 2), 
and tf^ = <5(±(e m _ 2 - £ m ) - ^ q c))- 



With the approximation gi — A ^ g 2 and the thermal 
average (n q ) = l/(e /3riw ' 3 — 1) one obtains 

W -2 = 77. "1C777 1 / 7^ d ±- 



4 ^ NMw„ 



± (e*^ - 1) 



(A8) 



As a next step the sum is replaced by an integral 
((l/A0E q ^ (a 3 /(27r) 3 ) / d 3 <j] and the density p = 
M/a 3 is inserted, 



w = A 2 s_ 2 f fq (q x - q y ) 2 + ±(q x + q y ) 2 



32tt 2 p 



(A9) 



After changing to spherical coordinates one gets 



W-2 



17A 2 .s_2 f 00 dq q 4 
1927rp J ~ q ± ( e ±/5^ - l) 



<J±. (A10) 



Assuming a linear dispersion relation cij g = cq, where c is 
the sound velocity, and using e = huj q = hcq one obtains 



W-2 



17A 2 s-2 



de 



1927rpc 5 ft 4 J "~ ± (ei/ 3 ^ 
= 17^ 2 S_2 (e m - 2 - £m) 3 
~ 192?rpc 5 ^ 4 e^( £ — 2 - £ ™) -l' 

In the same way we get 

17A 2 S +2 (£m+2 ~ £mf 



1) 



s± 



W- 



+ 2 



1927rpc 5 ft 4 e ^( £ -+2- £ -) - l' 



(All) 



(A12) 



with s+2 = (s — m)(s + m + l)(s — m — l)(s + m + 2). 

The transition rates for m — * m ± 1 can be calculated 
in the same manner as above with — 2 A « g%, 



W±i = 



A 2 s±i (e 



mil 



l2irpc 5 h A e ^( £ ™±i- £ -) - 1' 

l)(2m ± l) 2 , and p 



(A13) 



where s±i — (s ^p.m)(s ± m 
1.83 x 10 3 kg/m 3 E3 



APPENDIX B: LEVEL SPLITTING 

In this appendix we derive a formula for the tun- 
nel splitting energy which is applicable to potentials 
Krii,mi+i £ ^ with arbitrary Am = m t — m j+1 (m > 
to, > fii+i > m', i — 1, . . . , N — 1). According to Kato's 



theory^ the expansion of the resolvent 



G{z) 



1 



£f - AV 



(Bl) 



leads to a rigorous treatment of the perturbation theory, 
which is very useful to evaluate high-oijder perturbation 
terms. We use the notation of Messiah.ED 
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Let N be the order of the perturbation. Then the 
projection operator P = J2m l m ) ( TO I> consisting of the 
degenerate states {|m)}, and the operator (H — E^jP 
are expanded as follows: 

oo oo 

p = p o + E xNa(n) > ( h - E «) p = E xNb(n) , 

AT=1 



with 



A™ = 

B (N) = 



J2s kl VS k2 V---VS kN +\ 

(AT) 

E S kl VS k2 V ■■■VS kN+1 , 

(N-l) 



where 



S k 



-P , iffc = 0, 
%, if k> 1, 



1 



Qo — 1 — Ph — Qo- i '(ii- 

a fc ~ #o) 

and the sum 5Z(j\r) nas to be taken over all combinations 
fei, fc2, . . . , k]y+i with the restriction ki+k 2 + - • -+fc/v+i = 
AT. 

The following general secular equation must be solved: 



det (H a - X Ka) = det (C a ) = 0, 



(B2) 



where we have introduced the abbreviation C a = H a — 
xK a - The x are the eigenvalues of the perturbed states. 
H a and K a are defined by 



H a = PoHPPo = E° a K a + P E A^-BWPq, (B3) 

JV=1 

oo 

K a = P PP = Po + Po E ^AWPo. (B4) 

AT=1 

Thus we have now 

oo 

Ca = (El -x)Po + (E° a -x)E A^o^P) 

oo 

+ £ \ N P BWP . (B5) 



jv=i 



Equation (B5) is the general formula for finding the per- 
turbed eigenvalues and eigenstates. We apply it now to 
the situation of our two degenerate spin states \m) and 
\m'). The following derivation refers to the off-diagonal 
elements of Eq. ( p5|) . 

The factors PqW^Po and P P (Ar) P do not vanish if 
ki = k^+i = 0. As we look for the lowest-order pertur- 
bation that gives a contribution to the tunnel splitting 
E mm ', the projection operators S ki , i = 2,...,N, must 
not be equal to — Pq, i.e., ki ^ 0, i = 2, . . . ,N. Hence, 



we get the following combinations for the lowest-order 
perturbation, 

for A W : k 2 = k 3 = ■■■ = h- X = l,h = 2, 

h+i = ... = k N = 1, % = 2,---,N, (B6) 
for : k 2 = k 3 = ■ ■ ■ = k N = 1. (B7) 

In t he c ase of weak perturbation the second term of 
Eq. (B5) is much smaller than the third one. Thus the 
secular equation reads as follows: 

C a = (E° — x) Po + (diagonal elements) 



oo 



E AW E 



N=l 
N-l 

n £ 



v, 



(B8) 



Thus we arrive at formula ( p0|) 



APPENDIX C: APPLICATION OF 
KIRCHHOFF'S RULES 

In this appendix we make use of Kirchhoff 's rules (KI) 
and (K2) in order to evaluate the diagrams of the relax- 
ation paths. Each diagram and its evaluation is valid 
for the interval between two main peaks. The solu- 
tion t* n2 of the Kirchhoff equations between the peaks 
rii = m\ + m! x and n 2 = m 2 + m' 2 is not written down 
explicitly, since it is too lengthy and the calculation is 
straightforward. 

W^{A+im<H z <^(A + MB): 
From J = Ji + J 2 , Ji = J 3 + J 4 , Ji = J 5 + Js, -h + 
J 5 = Je,J3 + = Ji,Ji + Js = J, and AA" = 

J\T\ + J3T3 + J7T7, J3T3 = J4T4 + J5T5 + J 6 r 6 , J 2 t 2 = 

J\T\ + + J5T5, J 8 r 8 = J5T5 + Jqtq + J7T7, one can im- 
mediately evaluate the relaxation time r* 2 (H) = AN/ J 
(see Figs. |l| and 0). 
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FIG. 16. Spin relaxation paths (from m = 10 to m = —10) 
for J- (A + 13B) < H z < -2- {A + 345). Full lines: ther- 
mal transitions due to phonons. Dashed lines: dominant tun- 
neling transitions due to Bi and Hx terms. Dotted lines: 
tunneling transitions that lead to satellite peaks [included in 
the numerical diagonalization of the master equation (pq)]. 




FIG. 17. Serially reduced diagram associated with Fig. 
In order to understand the analytical evaluation of the re- 
laxation diagram in Fig. |l^ better, tunneling transitions that 
lead to satellite peaks are excluded. The relaxation times t„ 
are given in Eq. (p3[). 

From J = Ji+ J 2 , Ji = J3+J4+J7, J2+J4 = J5, J3+J5 = 
J 6 , J 6 + J 7 = J, and AN = Jin + J 3 t 3 + J 6 r 6 , J 3 r 3 = 

J4T4 + Jb T biJ2T~2 — J\T\ + J4T4, J7T7 = J4T4 + J5T5 + 

J6Tg, one can immediately evaluate the relaxation time 
t 2 * 3 (H) = AN /J (see Figs, pi and pi) . 




o> 



-2> 

-4) 
-6) 
-8) 




FIG. 19. Serially reduced diagram associated with Fig. [ii]. 
In order to understand the analytical evaluation of the re- 
laxation diagram in Fig. |l^ better, tunneling transitions that 
lead to satellite peaks are excluded. The relaxation times r„ 
are given in Eq. (E^j. 



(3) -2- (A 
From J = 



17B)<H Z <^(A + 40B): 

J1+J2, Ji — J3+J4, J4 = J5+J71 J2+J3+J5 — 
Je, Jq + J7 = J, and AiV = J\T\ + J 3 r 3 + J 6 r 6 , J 3 r 3 = 
J4T4 + J5T5, J 2 r 2 = J\T\ + J4T4 + J5T5, J 7 r 7 = J 5 r 5 + 
JqTq, one can immediately evaluate the relaxation time 
T 3 * 4 (Jf ) = AN/ J (see Figs. || and 




FIG. 20. Spin relaxation paths (from m = 10 to m = —10) 
for ^ (4 + 17B) < H z < (A + 40B). Full lines: ther- 
mal transitions due to phonons. Dashed lines: dominant tun- 
neling transitions due to B4 and H x terms. Dotted lines: 
tunneling transitions that lead to satellite peaks [included in 
the numerical diagonalization of the master equation 



FIG. 18. Spin relaxation paths (from m = 10 to m = —10) 
for -2- (A + 34B) < H z < -2- (A + 175). Full lines: ther- 

9MB v ' — — 9MB v ' 

mal transitions due to phonons. Dashed lines: dominant tun- 
neling transitions due to B4 and H x terms. Dotted lines: 
tunneling transitions that lead to satellite peaks [included in 
the numerical diagonalization of the master equation ([§8|)]. 



1G 




FIG. 21. Serially reduced diagram associated with Fig. |2C|. 
In order to understand the analytical evaluation of the re- 
laxation diagram in Fig. ^ better, tunneling transitions that 
lead to satellite peaks are excluded. The relaxation times r n 
are given in Eq. (E^j. 



APPENDIX D: FIRST-ORDER VS. 
SECOND-ORDER TRANSITION 

We show in this section that second-order transitions 
lead to a much faster relaxation of the spin system than 
first-order transitions if the coupling constants are equal. 
The relaxation rate I^ 1 ) of the cascade with transitions 
Am = ±1 has been calculated by Villain et alSJl (see 
Fig. H), 



3 l^i.ol 2 / \3 < 

2^7^ (e °~ £l) T 



3 IV, 



-PA 



1,0 



2tt h 4 pc 5 



- e -/3(eo-Ei) 
0A 



-/3A/s 2 



(Dl) 



A = 100 A is the energy barrier. 




|-9) 



FIG. 22. Cascade with Am 



-1 and H = 0. 



We have extended this expression by taking higher- 
order transitions into account. If we take a cascade with 
transitions Am = ±2, for the case s=10, we obtain 



r(2) = ±^H (£o - £2 ) 3 . 



-I3A 



2tt H 4 pc 5 

3 1^2,0 1 2 
2tt h 4 pc 5 



A 



1 3 



{s/2f 



j— /3(eo-E2) 
c-PA 



(D2) 



Comparing to the relaxation rate with s 
increase by a factor 



10, 



r(2) 



io 6 
"58 



= 64 



(D3) 



is obtained, assuming w (see Abragam and 

Bleaney,E3 p. 563, for experimental evidence). 

Now we calculate the relaxation rate by means of for- 
mula d|) with r« = 1/t*, Am = ±1, and = 1/r*, 
Am = ±2. If there is a fast transition via tunneling be- 
tween levels m = 4 and m' = —4 for H z = at T = 1.9 
K, we get the following more accurate estimation: 



rt 2 ) 
f(i) 



11.7. 



(D4) 



The same can be done if the fastest transition takes place 
via tunneling between levels m = 2 and m' = — 2 for 
H z = at T = 1.9 K, 



rt 2 ) 



49.0. 



(D5) 



From these results it is obvious that second-order transi- 
tions lead to a faster relaxation. Note that it is Eq. dl£ ) 
together with Eq. ( p^ ) which imply that the ratios (D4) 
and ( |P5| ) are of the same order as the ratio ( [D3| ) . This 
provides a theoretical justification for the approximation 
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